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1.0 INTRODUCTION 


This report summarizes the work conducted under this contract 
during the period of January 5, 1984, to April 4, 1984. Contract goals 
are to develop an efficient, comprehensive Si solar cell modeling program 
that has the capability of simulation accuracy of 5 percent or less. 

Computer simulation programs, in general, may be subdivided into a 
number of major tasks: (1) analytical method used to represent the 

physical system; (2) phenomena submodels that comprise the simulation of 
the system; (3) coding of the analysis and the phenomena submodels; 

(4) coding scheme that results in efficient use of the CPU so that CPU 
costs are low; and (5) modularized simulation program with respect to 
structures that may be analyzed, addition and/or modification of phenomena 
submodels as new experimental data become available, and the addition of 
other photovoltaic materials. These tasks are briefly discussed below. 

Computer modeling simulations have been shown to be very useful in 
the development of semiconductor devices in those cases where the simulation 
is an accurate representation of the physical device. However, to be an 
effective aid to the experimentalist and to become an equal partner in 
the technologies used in device development, it may be required to 
operate the computer program frequently each day in an active develop- 
mental program. For frequent use, as required in solar cell development, 

CPU costs must be low. Moreover, low CPU cost allows for engaging in 
computer experiments, which can be made to be a very useful and powerful 
technique. 

Computer modeling using numerical integration methods in Si device 
technology have usually shown fair-to-good agreement with experimental 
data. However, CPU costs for the execution of computer programs that 
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are based on numerical integration methods are typically high for their 
use as a laboratory or manufacturing tool [1]. The number of bias 
points that are required to study optimized device designs is usually 
more than 1,000 to 2,000 runs. Similarly, a comprehensive study involving 
device structures or new types of devices typically exceeds 2,000 runs. 

In most cases, the cost of such modest studies for the benefits gained 
is not attractive. It is clear that conducting in-depth studies is not 
cost-effective; for this reason, numerical integration methods are not 
used for comprehensive studies. 

A method has been developed that, in principle, should result in 
simulation accuracy equivalent to that demonstrated by numerical integration 
methods, but with low CPU costs. This method uses recursion relationships 
to solve for the constants of integration and is subsequently described. 

Simulation accuracy is determined by both the accuracy of the 
analytical method representing the device and the accuracy of the phenomena 
submodels in representing the corresponding experimental data related to 
material properties. For most efficient use of the CPU, the accuracy of 
the analytical method and of the phenomena submodels should be equivalent. 
For example, even if the analytical method accurately represents the 
device, simulation results may not agree with experimental data if the 
phenomena submodels are accurately represented. The reverse is also 
true. In solar cells, the phenomena submodels that produce first-order 
effects in terminal characteristics are: absorption curve, built-in and 

induced electric fields, bandgap, lifetime, mobilities, diffusi vities , 
photoexcited carrier concentrations, surface recombination velocities, 
junction transport, etc. Some of these will be discussed in this report. 

The representation of the phenomena submodels must take on an importance 
equal to the analytical method used to represent the system. 
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Codings of the analysis and phenomena submodels are separated in 
the program. Moreover, the coding of the submodels is modularized so 
that there is maximum flexibility to add, subtract, and/or modify them. 

The efficiency of the coding may only be demonstrated by the accuracy 
of the simulation results with experimental data and the corresponding 
CPU execution time. 

The spectrum of structures and cases that may be treated by the 
simulation program depends on the analytical method or algorithm used 
to represent the device, generality of the boundary conditions imposed, 
comprehensiveness of the phenomena submodels, and the generality of 
representation of each of the submodels. These issues will be described 
below. 
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2.0 DESCRIPTION OF THE INVESTIGATIONS CONDUCTED 

2.1 Analytical Method Used to Represent Solar Cell Structures 

RT1 has been developing computer modeling programs for semiconductor 
devices since 1976 [2-9]. The major element that has evolved out of 
these modeling programs is the analytical method that is used to analyze 
physical systems. To obtain an accurate representation of a physical 
system, the system geometry is divided into a number of segments (mesh 
points). In a one-dimensional geometry, the segments are defined by a 
series of parallel planes. The separation of the planes defining the 
segments determines the simulation accuracy. The segments may be taken 
as thin as required to obtain the accuracy desired. The transport 
equations, governing the behavior of the physical system, are applied to 
each segment, and a closed-form solution is obtained in each of the 
segments. It has been demonstrated that when applying the boundary 
conditions to the solutions at each of the planes, there exists a recursion 
relationship between the constants of integration of the solutions 
obtained in each segment. 

This method has been applied to semiconductor devices and has 
been shown to give excellent agreement with experimental results. It 
has been applied to devices for which the solutions are functions of 
position and time, such as photodiodes [8] and avalanche photodiodes 
(APD) [8], The method has also been applied to solar cells, where the 
solutions are functions only of position [unpublished]. 

Convergence problems do not impose any conditions on the mesh 
point separation. The segments may be made arbitrarily thin or 
arbitrarily large. For maximum accuracy in simulating a physical 
system, a large number of segments may be required. The segments may be 
increased beyond 1,000 for greater accuracy, however, it has been found 
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that this is not always necessary. In the case of APDs, a device for 
which accuracy of representation is difficult, simulation results were 
studied for f = 2 up to f = 80 to determine the agreement between modeling 
and experimental data. It was found that the results obtained for f = 

15 were within 1 percent of the results obtained for f = 80 [8]. This 
clearly shows these real advantages of the recursion relationship method 
over the numberical integration method: (1) a "guess solution" is not 

required; (2) convergence is always obtained for any number of mesh 
points; (3) a smaller number of mesh points may be used; (4) transport 
equations are solved by means of recursion relationships; and (5) CPU 
running time is reduced significantly. 

During this contractual period, the analytical method described 
above has been applied to silicon solar cell structures; however, it has 
been generalized so that it may be applied to heterostructures encountered 
in amorphous silicon, III-V and 1 1 -VI material -based cells. Cascade 
cell structures may also be analyzed by means of this solution. The 
silicon solar cell structures that may be analyzed are given in Table 1, 
including the effect of an oxide-charged insulator (OCI). 

Final results of the recursion relationships have not yet been 
obtained and will be reported in the second quarterly report. 

2.2 Analytical Representations of Phenomena Submodels 

In this section, a brief discussion is presented of some of 
the phenomena submodels and the representations used in the simulation 
program. 


2.2.1 Mobility 

Carrier mobilities in silicon have been studied by Dorkel and 
Leturcq [10], denoted DL mobility. They have assessed experimental and 
theoretical studies and developed analytical relationships which represent 
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Table 1 


Silicon structures that may be studied 
using. the simulation program. 


6 







low field electron and hole mobility, depending on impurity concentration 
(ionized impurity scattering), temperature (lattice scattering), and 
injection level (carrier-carrier scattering). The DL mobility formulation 
has been selected as the primary method to be used in the simulation 
program, because the mobility components that arise from the three 
scattering phenomena are easily discerned. In addition, the composite 
form of the DL mobility relationships are relatively simple to apply, 
and it contains the mobility component arising from carrier-carrier 
scattering which becomes important in degenerate material and/or high 
injection levels. However, it appears that the DL relationships may not 
accurately represent mobility data for net impurity concentrations above 

I O O 

2-3 x 10 cm under low injection levels and for nonequilibrium values 

36 *“6 

of the electron-hole product (p n or p n ) above 10 cm . This 

r r n n r p p 

aspect will be dealt with subsequently. 

The DL mobility relationships, which provide for silicon 
mobility dependency on impurity concentration, temperature and injection 
level are given by the expressions: 
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1.025 


1 + ^ 1 . 68 ^ 


1.43 


- 0.025 


(5) 


In the above set of equations, both electron and hole mobilities are 
represented and for which the baseline constants are listed in Table 2. 
Using these relationships, both electron and hole mobilities may be 
obtained in n-type material where N represents the net donor concen- 
tration, Similarly, electron and hole mobility may be determined 

in p-type using the appropriate mobility parameters in Table 2. 

Eqs. (1), (2), and (3) represent the mobility components 
produced by lattice, ionized impurities, and carrier-carrier scattering. 

The latter is important in degenerate material and at high injection 
levels, and is usually negligible for most other cases. On the other 
hand, the lattice scattering and ionized impurity scattering are determined 
under conditions close to thermal equilibrium. 

Lattice scattering mobility, Eq. (1), is produced by a scattering 
phenomenon attributed to accoustic phonons. This component decreases 
with increasing temperature, and the parameters, y^g and "a," depend on 
the carrier type and the temperature range. It is this component which, 
in large measure, determines the temperature behavior of mobility. It 
is clear from Table 2 that y^g > U|_pcr 

Ionized impurity scattering mobility, Eq. (2), is produced 
because of scattering of carriers by immobile ionized impurities. This 
mobility decreases with increasing impurity concentration or decreasing 
temperature. Implicit in this formulation is that the effect of scattering 
is the same for electrons and holes in the presence of ionized donors as 
well as for ionized acceptors. 
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Parameter 

Electrons 

Holes 

m Lo 

1430 cm 2 (V sec) -1 

495 

a 

2.2 

2.2 

A 

4.61 x 10 17 (cm V sec K 3/2 ) -1 

1 x 10 17 

B 

1.52 x 10 15 cm -3 K -2 

6.25 x 10 14 


Table 2. Baseline values of mobility parameters using the 
Dorkel-Leturcq [10] representation. 
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Carrier-carrier scattering mobility, Eq. (3), becomes important 

where the electron-hole product is greater than 10 ^ cm ® for 10^ cm ^ 

32 -6 

impurity concentration and when the product is greater than 10 cm 
19 -3 

for 10 cm impurity concentration. These conditions occur for injection 

levels produced under solar concentrations exceeding 10 to 20 suns. 

Electron mobilities predict experimental values to within 

5 percent in the temperature range above 200 K, and for doping levels 
18 - 3 

less than 3 x 10 cm" . At 300 K, hole mobilities are predicted to 

1 8 -3 

within 5 percent or less, and for doping concentration less than 2 x 10 cm . 

There is agreement with experimental data in the temperature range 200 

to 300 K. However, hole mobility data above 300 K and the DL mobility 

relationships cannot be assessed at this time [10]. 

Due to limitations of applicability of the DL formulation in 

some temperature and doping concentration ranges, a second representation 

has been included in the mobility submodel. This representation is due 

to Arora et al. [11]. It is denoted as the AHR mobility formulation. 

Composite mobility relationships for electrons and holes have 

been obtained empirically, as a function of temperature (lattice scattering 

mobility) and impurity concentration (ionized impurity scattering mobility), 

based on experimental data and the Brooks-Herring theory of mobility. 

The relationships predict electron and hole mobility to within ±13 percent 
20 -3 

up to 10 cm doping concentration and in the 250 to 500 K temperature 
range. The relationships are given by: 
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The above formulation may be used in those temperature and doping concen- 
tration values where the DL formulation exhibits a poor representation 
of experimental data. However, Eqs. (6) and (7) only apply to low 
injection level cases, in which case we will resort to the DL formulation. 

The DL [10] and AHR [11] relationships are compared in the 

■14 19 -3 

impurity concentration range of 10 to 10 cm in Figures 1, 2, and 3 
for 200, 350, and 500 K, respectively , under low injection level, 
pn = 10 cm - . Agreement for electron mobility between the formulations 
is excellent at 350 and 500 K and not as good at 200 K. The overall 
agreement for holes is not as good as it is for electrons. However, 
these results show that either formulation may be used to represent 
mobility data under low injection levels. 

2.2.2 Diffusivities 

In most cases, the carrier diffusivities used are of the 
Einstein form represented for electrons by 




( 8 ) 


However, in degenerate and/or high injection levels, the Einstein 
relationship is a poor approximation. 

Consider a degenerate p-type semiconductor under thermal 
equilibrium (i.e., open-circuit conditions) as shown in Figure 4(a). 
Under these conditions, the electron current, given by 
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200 K 



12 


Figure 1. Comparison between the DL and AHR relationships for electron and hole mobility 

at 200 K under low injection level. 
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Figure 2. Comparison between the DL and AHR relationships for electron and hole mobility 

at 350 K under low injection level. 


500 K 
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Doping Concentration, cm 

Figure 3. Comparison between the DL and AHR relationships for electron and hole mobility 

at 500 K under low injection level. 
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is zero, from which the electron diffusivity is given by 


D - "n n P° E x 
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dn 


(10) 


£0 


dx 


E x in this case is the built-in potential that is established to produce 
an equal, but opposite drift current to the diffusion current. The 
built-in field is obtained from a relationship for the conduction bandedge, 
E c . The bandedge expression may be placed in a form given by 


E = -i)> + constant 
c 


(ID 


where ip is the electrostatic potential. The built-in potential is 
obtained from Eq. (11), represented by 


: .5 

'x dx 


d\p 

dx 


( 12 ) 


To obtain dn nQ /dx, we use the relationship 


dn ™ dE„ 

po _ po c 

dx “ aE„ dx 
c 


(13) 


and substituting from Eq. (12), the result is 


dn an„„ 

M = e — 2° 


dx 


"x 3E. 


(14) 
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The relationship for electron concentration (using the rigid band 
approximation) given by 


n 


po 


N c F l/2^ n c 


( 15 ) 


where N c is the density-of-states at the conduction bandedge, F-j ^ is 
the Fermi -Dirac integral of 1 /2-order, and 


n Q = q 


E Fo‘ E c 

kT 


(16) 


is the argument for thermal equilibrium. Using the relationship 
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^ F l/2 (n o^ 


it can be shown then [12] that 


(17) 


an 0 [ F l/;2 (n o^ = F -l/2^ n o^ ’ 

where F_-|/2^ n o^ 1S Fermi-Dirac integral of the -1/2-order. Substituting 
Eqs. (14), (15), (16), (17), and (18) into Eq. (10) results in the 
electron diffusivity in degenerate material given by 
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i/2 ( y 

-l/2 (n o ) 


(19) 


A similar procedure and relationship may be obtained for the case of 
high injection level, by elevating temperature or by another method in 
which the same electron-hole pair concentration is produced as through 
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irradiation by a high photon flux as shown in Figure 4(b) under dynamic 
equilibrium and open-circuit conditions. The latter assumes that a 
uniform distribution of electron-hole pairs is generated. The diffusivity 
is given by 


D 


n 



F 

F 


1/2^V 

-l/2 (n n^ 


( 20 ) 


for which 



( 21 ) 


where Ep n is the quasi-Fermi level for electrons. Under thermal equilibrium 

conditions, Eq. (20) reduces to the relationship given in Eq. (19). 

Corresponding relationships exist for holes. 

Figure 5 shows the electron and hole diffusivities, normalized 

with respect to the Einstein diffusivity, plotted as a function of their 

respective arguments (n and -eg-n)' which contain the Fermi level. 

Increasingly positive arguments represent increasing degeneracy for 

which the diffusivities show significant increases. Negative arguments 

represent nondegenerate cases where the diffusivity ratios approach 

unity. For argument values of +10, the diffusivity ratios are greater 

than six. This represents a condition where the Fermi level penetrates 

20 -3 

the corresponding bands by 0.15 eV, which are attained for 10 cm 
doping under low injection. High injection adds to this penetration, 
and the diffusivity ratio increases still further. 

2.2.3 Boundary Conditions at Depletion Region Edges 
Boundary conditions imposed at depletion region edges are 
briefly discussed below. Due to the controversy that may surface when 
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I J I 



Figure 4. Physical systems used to describe carrier diffusivities: (a) degenerate 
semiconductor in thermal equilibrium; and (b) degenerate and/or high 
injection level in a semiconductor in dynamic equilibrium. 
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Figure 5. Electron and hole diffusivities normalized with respect to the Einstein 
relationship for degenerate and/or high injection level conditions for all temperatures. 
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discussing these boundary conditions under high injection levels, special 

attention was devoted to this subject [13-19], During the course of 

this investigation, a number of significant results were obtained, and 

these will also be presented below. 

Consider the geometry where Xg and x^ denote the depletion 

region edges in the n- and p-regions, respectively. The symbols for the 

corresponding carrier concentrations are denoted P n (x 2 ^ 5 n n^ x 2^’ p p^ x 3^’ 

and n (xo). In general, the relationships that exist between these 
P ^ 

concentrations are given by 

_!!bi 

P n ( x 2 ) = P p ( x 3^ e kT ’ (22) 

for holes, and for electrons by 

n p ( x 3 ) = n n ( x 2 ) kT » (23) 

where V... is the nonequilibrium built-in potential established in the 
depletion region. While the expressions in Eqs. (22) and (23) are valid 
for both equilibrium and nonequilibrium conditions, at high injection 
levels the quantities P p ( x 3 ) an ^ n n ( x 2 ) are unknown. This gives rise to 
two different forms of the boundary conditions for low and high injection 
levels. The different forms are a direct result of the charge neutrality 
conditions that exist under low level injection, given by 

n n^ x 2^ n no^ x 2^ = V x 2^ (24) 

in the n-region, and in the p-region by 
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p p< x 3> Ppo 1 ^ 1 = N A (x 3> • 


(25) 


where Nq(x 2 ) and N^(x 2 ) represent donor and acceptor concentrations in 
the n- and p-type regions, respectively, at the depletion edges. 
Substituting Eq. (24) into Eq. (23), and Eq. (25) into Eq. (22), gives 
the boundary conditions for low injection level: 

p n (x 2 ) " Ppo 1 ^' 6 ^ ’ (26) 


and 


n p (x 3> 


n no( x 2> e 


pV bi 

kT 


(27) 


The nonequilibrium built-in potential is defined by 


V bi ■ V bio - V J 


(28) 


where V^. 0 is the built-in potential at thermal equilibrium, and Vj 
is the photovoltage. Therefore, under low injection levels, p n (x 2 ) an ^ 
np( x s) may be determined. 

In the high injection case, the charge neutrality condition 
imposed in the n-region is given by 


n n (x 2 ) - P n (x 2 ) = N d (x 2 ) 


(29) 


and in the p-region by 
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where P n ( x 2 ^ ma ^ ^e com P 3r3 b1e to or even greater than Np^) and 
similarly for n p (x 3 ) and N A (x 3 ). One form of the boundary condition 
under high injection was first reported by Fletcher [14] and given by 


p„(x 2 ) = 




1 - t 


(31) 


for holes at x 2 and by 


n P <*3> = 


[N d (x 2 ) + N a (x 3 ) 5 ] K 

T? 


(32) 


at x 3> where 

qVj 

= n i e ( x 2’ T,i y n ie( x 3 ,T ’V FT 
V x 3^ V x 2^ 


(33a) 


The quantity Vj is the photovoltage generated across the depletion 
region, x 3 -x 2 - The form of f; is somewhat different from the results 
obtained in the literature because it includes the effects of bandgap 
narrowing through the intrinsic concentrations n i - e (x 3 ,T,Np) and n -j e ( x 3 »T >N^) . 
To show this more explicitly, Eq. (33a) may be placed in the following 
form: 


P [AE r (xJ + AEp(x,)]/2kT 
n^ 0 (T)e 62 G 3 

N A^ X 3^ N D^ X 2^ 



(33b) 


This form shows that K is dependent on the sum of the bandgap narrowing 
of the depletion edges. If the doping concentration at x^ and x 3 are 
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1 x 10 cm , the sum of the bandgap narrowing is 0.086 eV and e is 
increased by a multiplicative factor of 28 compared to pure material. 

An alternative form of these boundary conditions is attributed 
to Misawa [13]. While the Fletcher formulation is based on quantities 
defined at the depletion region edges (x£ and x^), the Misawa formulation 
also includes quantities defined at the outer boundaries of the n- and 
p-regions, which represent the ohmic contact interfaces (x = 0 and x^) 
in a solar cell structure. The two relationships are given by 


P n (x 2 ) 


Ppo (x 5> e 


-9 T [V(x 2 )-V(x 5 ) + |; 


qu P P P 


dx] 


(34) 


and 


= "no (0)e 


-3 T [V(0)-V(x 3 ) + }^A-dx] 


(35) 


where p pQ (xg) and n no (0) are the acceptor and donor concentrations at 
the ohmic contacts, and V(x 5 ) and V(0) are the corresponding potentials. 
The disadvantage of the Misawa form is the photocurrent J p and J p are 
required to be known in order to calculate p n (x£) and n p (xg). 


2.2.4 Boundary Conditions Imposed at the Mesh Points 
The analytical method described in Section 2.1 solves, in closed 
form, for the minority carrier concentrations at each mesh point. Therefore, 
each of the closed form solutions contains two constants of integration, 
which solutions must be obtained through the imposition of boundary 
conditions of each of the mesh points. 
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are: 


In the n-region', the boundary conditions at the j 


th 


mesh point 


P . (y . ) = p 

K nej'- , oi / K 


ne( j+1 ) 


(AE ./kT) 

(0)e VJ 


(36) 


and 


J pj^ y oj^ J p(j+1) 


( 0 ) 


(37) 


where P ne j(y Q j) and p n e(j+l)^ represent the photoexcited hole concen- 
tration, J pj (y oi ) and ^p(j+i)(0) the h°^ e current density, y Q ^ the 

X L x L 

separation between the j -1 and j n mesh point, and aE . the discontinuity 

’ J 

in the valence bandedge. If AE yj . > 0, a discontinuity in hole concentration 
exists and provides for carrier confinement. A corresponding set of 
conditions are imposed in the p-region. 

2.2.5 Bandgap Narrowing and Intrinsic Concentration 
Bandgap reduction has been shown to occur in silicon with 
increasing impurity concentration, which is attributed to the broadening 
of impurity levels and which ultimately overlap with the conduction and 
bandedges when donors and acceptors, respectively, are involved. This 
effect is well established from the results of band structure analysis 
[20-22], and has been shown by direct means through the measurement of 
the intrinsic optical absorption threshold [23], The largest narrowing 
observed is 0.068 eV, using this technique, in n-type silicon for a 
donor concentration of 9 x 10 cm [23]. These experimental results 
were employed to develop models by a number of investigators to obtain 
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improved agreement between experimental data and theoretical calculations 
related to effects produced by high impurity concentrations in the 
emitter and base regions of bipolar transistors [24-26]. 

More recently, other methods have been used to measure bandgap 
narrowing [27,28]. These data have been analyzed and empirical relation- 
ships obtained which describe the bandgap and bandgap narrowing with 
increased impurity concentration and which are denoted here as the SD 

[27] and D [29] relationships. In most cases, the experimentally deter- 
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mined bandgap narrowing in the 10 to 2 x 10 cm range is greater 
than the values resulting theoretical analysis [30-32]. The Slotboom 
and de Graaff (SD) values have been used to calculate injected electron 
current under contacts in logic circuits [33,34], to obtain improved 
agreement in calculating built-in potentials [35], and to model solar 
cells [36]. The results appear to be more encouraging for the SD than 
for the D relationships; therefore, in the following discussion the SD 
results are used. 

Slotboom and de Graaff measure the bandgap by two indirect 

methods that involve the I^-V^g characteristic in NPN bipolar transistors, 
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as a function of impurity concentration (4 x 10 to 2.5 x 10 cm ) 
and temperature (150 K to 400 K). Following MacFarlane et al . [37], a 
linear approximation is assumed for the bandgap above 250 K and represented 
by 


E g (T,N) 


'GO 


aT - AE g (N) 


(38a) 


where a = 3.855 x 10 ^ eV(K) ^ is the temperature coefficient, E 


GO 


1.206 eV 


is a constant, and AEg(N) is the bandgap narrowing. Both E^q and a are 

independent of impurity concentration and temperature, while aE~(N) is 

b 
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dependent only on the net impurity concentration. If the bandgap, 
E~(T,N), is measured, the bandgap narrowing may be determined from the 
relationship 


AE G (N) = E qo - 3.855 x 10" 4 T - E fi (T,N) 


(38b) 


In analogy with nondegenerate material, Slotboom and de Graaff use the 
relationship 

o -E r (T ,N)/kT 

n ie (T ’ N) = N c N v e (39) 


2 

to calculate Eg (I , N ) once they obtain n^ (T,N) from the measurement of 
collector current. N c and N y are the density of states at the conduc- 
tion and valence bandedges, respectively, and are assumed to be invariant 
in form and in their values from those in pure silicon. These assumptions 
constitute what is commonly referred to as the Rigid Band approximation, 
because the assumptions are equivalent to assuming that the bandedges 
remain parabolic even in degenerate cases. Although this is almost 
certainly not true, it still represents an approximation that is more 
realistic than the estimates that are required in the bandstructure 
calculations [20-22]. As a result, Eq. (39) may be written in the form 


n- e (T,N) 



AE r /kT 

(T)e G 


(40) 


where the square of the intrinsic concentration in nondegenerate material 
is given by 


n? 0 (T) = 9.6 x 10 32 T 3 e Eg ° AT 


(41) 
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Slotboom and de Graaff state that their data demonstrates that AEg(N) 
is independent of temperature. Therefore, they measured n^ (T,N) only 
at 300 K over a range of impurity concentration values. It is then 
assumed that the form of Eq. (40) applies to the doping and temperature 
studied. An empiracle relationship is obtained to represent bandgap 
narrowing and is given by 


aE g (N) = 9 x 10 

Substituting for net concentration values, N gives the bandgap narrowing 
obtained from the measurement of n^ g (T,N). 

The Slotboom-de Graaff assumptions may be questioned. Although 
they have maintained electron injection into the degenerate base region 
at a low level, they claim that this is reason enough to represent the 
pn-product by Boltzmann statistics. The temperature behavior of hole 
mobility in the base region is determined by measuring the sheet resistance 
of the base region. This hole mobility temperature dependency is then 
ascribed for electrons as well. Moreoever, the temperature dependency 
of bandgap narrowing is inferred to be constant as a result of measuring 
the temperature behavior of the emitter-base voltage (V^g) and of the 
electron mobility. The validity of the above assumptions is not questioned 
in general, but only as it applies to determining small effects such as 
bandgap narrowing. 

A more realistic method to determine aEq(N) from measured 
values of n 1 - e (T,N) is to use Fermi-Dirac statistics when studying effects 
in degenerate materials. In the following, we present the results of a 
study that applies Fermi-Dirac statistics to bandgap narrowing. In 
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place of Eqs. (39)-(42), we substitute the relationship that applies to 
degenerate and nondegenerage cases given by [12]: 


.2 /T _ 2 F l/2 (n)F l/2 ( ~ c G~ n) 


n- Q (T,N) = n‘ 


le 


10 F l/2 tn i ,F l/2 ( " e G“V 


(43) 


where F-j^x) represents the Fermi-Dirac integrals of the 1/2-order, and 


e f- e c 
n = IT 


(44) 


E Fi" E C 

kT 


E c- E v 


£ G = kT 


(45) 

(46) 


Ep, E Fi , E c , and E y are used to represent the usual parameters. 

The bandgap value, Eg(T,N), is determined by substituting the 
Slotboom-de Graaff experimentally determined values of n i - e (T,N), and the 
calculated values of n io (T) and n- into Eq. (43) through a computer 
subroutine that requires that the relationship be satisfied. Eg(T,N) 
is determined by this method over 1 x 10^ to 3 x 10^ cnT^ concentration 
range and 250 to 500 K. temperature range, where the invariance of aE^(N) 
with temperature is also assumed. The bandgap narrowing obtained from 
the two methods (SD method and the method used in this study) using 
Eq. (38b) are shown in Figure 6. These results show that for concentration 
values less than 1 x 10 cm - , both methods give identical results over 
the temperature range. However, the method based on Fermi-Dirac statistics 
(F-D) results in larger bandgap narrowing values compared to the results 
based on Boltzmann statistics (B), and the difference increases with 
increasing concentration values. Even though the experimental values of 


28 


Bandgap Narrowing, A E Q , meV 


Legend 

Fermi-Dirac Statistics (F-D) 

Boltzmann Statistics (B) 




» Figure 6. Bandgap narrowing calculated using Fermi-Dirac and Boltzmann 
statistical forms, and based upon nj„ (T,N) calculations using the Slotboom and 

deGraaf data at 300 K. 


29 


n i e ( T ’ N ) obtained by Slotboom and de Graaff for which aEq(N) is independent 

of temperature, the F-D method nevertheless shows that AEg(N) does 

exhibit a temperature dependency. Figure 7 shows the intrinsic concentration, 

n. (T,N), versus 1/T for three concentration values. It is seen that 

19 -3 

n^ e (T,N) is higher by approximately one decade for 3 x 10 cm compared 
to 1 x 10 17 cm 3. 


2.2.6 Carrier Lifetime 

It is well known that carrier lifetime plays a major role in 
obtaining high efficiency in silicon solar cells. As a result, a 
number of studies are contemplated that involve effects produced by the 
lifetime behavior which have not yet been resolved. For these reasons, 
composite lifetime relationships have been incorporated as submodels in 
the computer simulation program. The composite relationships contain 
recombination terms representing the Shockley-Read-Hal 1 , trap-assisted 
Auger, and band-to-band Auger processes. 


2. 2. 6.1 Shockley-Read-Hal 1 (SRH) 

Recombination processes that apply at low injection levels are 
referred to as the Shockley-Read-Hal 1 (SRH) recombination [40,41]. This 
recombination relationship may also be applied to Auger recombination 
through trap-assisted processes. In silicon, the latter is most important, 
because silicon is characterized by indirect optical transitions. 

The SRH recombination-generation relationship, which involves 
the electron and hole capture and emission at low injection levels through 
a single trap level, is given by [40,41]: 


Pn " 


n 


ie 


SRH (E.-Ep . )/kT 

T p[ n+n ie e 


r n [ P +n if 


■( E t- E Fi> 7kT - 


(47) 
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Intrinsic Carrier Concentration, nj e (T, N) 




The minority carrier lifetimes (t and t ), carrier concentrations (p 

j p n 

and n), intrinsic concentration (n ie ), trap energy level (E t ), and the 
intrinsic Fermi level (Ep^) may all be sensitive to position, temperature, 
impurity concentration, and, for a number of these parameters, to injection 
level . 

Kendall [38] has shown that experimental hole lifetime in as- 
grown n-type silicon ingots may be described by the expression - 


(Nn> = 


V"d 




1 + 


Do 


(48) 


where Np(x) is the donor concentration which may be position-dependent. 

The parameter Xp Q is the lifetime in material for which Np << V and 

N n is the threshold of donor concentration where recombination concentration 
Do 

begins to increase significantly. These parameters are listed in Table 3. 

By analogy, Fossum [39] has assumed a similar relationship for electrons 
given by 


T n ( V 


no 


1 + 


Ao 


(49) 


where x and N. have similar meanings as in Eq. (48). In those regions 
no mo 

for which Np and are position-sensitive, the minority carrier lifetime 
will also be position-dependent. 

Eqs; (48) and (49) may be made more general where the form of 
the relationships of Xp Q and x nQ are given by 


X 


0 


1 


crv 


th 



(50) 
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Table 3. Summary of lifetime constants. 
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The parameters a and are the recombination cross section and the 
thermal velocity, respectively, and N tQ is the recombination concentration 
when N a = N Aq or N Q = N Qo in the expression [42]: 


N t (x) = 1 


to + N(x) 


(51) 


In addition, the recombination cross section is represented as a function 
of temperature by [42]: 


a 


a(300K) ( 


300 

T ' 


(52) 


Using the constants provided for x , a(300K), and v^ in Table 3, N^ Q 
may be calculated and is also listed. Minority carrier lifetimes represented 
in Eqs. (48) and (49) are now expressed as functions of impurity concen- 
tration as well as temperature. 

2.2. 6. 2 Trap-Assisted Auger (TAA) 

Trap-assisted Auger (TAA) [43-47] is distinguished from the 
band-to-band Auger recombination, because TAA recombination is sensitive 
to trapping level concentration. As a result, TAA may be sensitive to 
fabrication processes. This arises in those materials and/or from 
fabrication processes that use high temperatures in which shallow trap 
concentration is influenced. The trap concentration model used for TAA 
is represented by Eq. (51), which was introduced for use with SRH 
recombination. 

The composite recombination model formulated considers that a 
single trap level contributes to both TAA and SRH recombination [42,48]. 

This is consistent with the single-level analysis that has been extensively 
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used in the literature and for which there is ample agreement between 
experimental data and simulation results. Although this analysis, 
represented by Eq. (47), was developed using the SRH capture and emission 
processes, it is reasonable to extend this relationship to include TAA 
recombination [42,48]. However, in SRH the lowest lifetime is obtained 
for those cases where the trap level ^ Ep^ [40,41], while in TAA it 
is obtained where E^ ^ E c in n-type and E^ ^ E y in p-type [42], 

The relationship for hole lifetime in n-type is given by 

[42,48]: 


x pTAA 


T p N tn"„ 


and for electrons in p-type by 


(53) 


x nTAA T N p 
n tp r p 


(54) 


where T and T are the TAA recombination coefficients for holes and 
p n 

electrons, respectively, and n n and Pp are the nonequilibrium majority 
carrier concentrations. Hole and electron recombination coefficients 
may be represented by the expression [43]: 


T = 2.23 x 10~ 26 E |5/2 

u/) 2 (E g (T,N)-E') 3/2 eJ!(T,N) 


where E' = E c ~E t for holes in n-type and E 1 = E^-Ey for electrons in p- 

★ 

type. The quantities e r and m are the relative dielectric constant and 
the appropriate effective mass. Eq. (55) predicts that electron recom- 
bination increases as E.-E w decreases, and is lowest for E.-E = E~(T,N). 

C V t V b 

A corresponding statement may be made for hole recombination. 
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2. 2. 6. 3 Band-to-Band Auger (BBA) 

A band-to-band Auger process may also be present in silicon, 
which does not require an intermediate trap level [49,50]. This recom- 
bination process takes place in degenerate material and/or under high 
injection levels. The recombination relationship is given by [42]: 


R BBA = (B n n + 8 p p)(np ‘ n ie> 


(56) 


where B and B are the recombination coefficients for electrons in p- 
n p 

type and holes in n-type, respectively. While the values given for the 
coefficients in the literature show a significant spread, those listed 
in Table 3 appear to be the most reasonable. 

2. 2. 6.4 Composite Relationship 

The composite relationship is obtained by combining the SRH 

and TAA recombination rate into Eq. (47) and adding Eq. (56) to the 

result. In Eq. (47) the hole and electron lifetime, x and x , are 
^ p n 

obtained in a form that includes SRH and TAA using the form [42,48]: 


1 

x 



J 

t TAA 


(57) 


for holes and electrons. This results in the relationship 


X 



LT p (x)n n (x) + ° p v thp JN tn tx) 


for holes in n-type using Eqs. (48) and (53), and 


(58) 


T n (x) = [T (x)p 


CT n v thn 




(59) 
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for electrons in p-type using Eqs. (49) and (54). After Eqs. (58) and 
(59) are substituted in Eq. (47), the result is represented by r srH/T AA' 


When Rg B ^ and r srh/TAA are a< ^^ ec * 5 the electron lifetime may be obtained 


given by 


Bn + Bp + 
n p p K p 


n - +n ie e 


TE^TTkT 


P 1 
" T n 

u r\ 


V n ie e 


-( E t- E Fi )/kT 


.. +a v.. jN. „ 

p p p thp J tp 


~T p +a V. . JN ~ 

• n r p n thn J tp 


p | , i + p e° + ^ 

K pe I p n 

H ' K pe pey 


where 


(60) 


and 






( 61 ) 

(62) 


The subscripts o and e represent equilibrium and excess carrier concen- 
trations, and n p and p p represent the total, nonequilibrium electron and 
hole concentrations, respectively. A similar relationship may be obtained 
for holes. 


2.2.7 Induced Surface Electric Field Due to PCI 
The simulation program provides for the option of imposing an 
immobile charge on the irradiated surface of the solar cell. This charge 
may be distributed arbitrarily throughout the insulator [53-63]. Never- 
theless, from Gauss' law the total charge may be considered to reside at 
the interface. An aiding field in an n-type surface, which reduces 
surface recombination and improves the junction collection efficiency, 
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requires the charge to be positive. This tends to repel minority carrier 
holes and to produce an electron accumulation layer. The model for this 
structure is shown in Figure 8. 

A number of methods have been used to obtain the relationships 
for the electric field intensity and the electron accumulation distribution 
[53-63]. The analysis follows from Poisson's equation for the one- 
dimensional case 

dx^ = e £ n d ( x ) + ' n n ^ x ^ ‘ ( 63 ) 

The boundary conditions imposed on the solution for E(x) are 

Q s 

E (0) = ~ (64) 

x e 


and 


E x (x 2 ) - 0 


(65) 


Integrating Eq. (63) results in the relationship for the electric field 


E s ( x) = | ] tn n (x) - p n (x) - N D (x)]dx 


( 66 ) 


This field component is added to the built-in electric field relationship 
that is discussed in Section 2.2.7. In Eq. (66) the mimobile surface 
charge density is given by 


j: 


[n n (x) - p n (x) - N p (x)]dx 


(67) 


and satisfies Gauss' theorem. 
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Depletion , 
Region 
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A second method may be used to obtain a relationship for the 
electric field, which is an approximation. The electric field is given 
by 


E 


s 




1 - e 


-y(* 2 ) 

e 

-y(x 2 ) 


( 68 ) 


where 


Y = [n n (0) - N Q ] , (69) 

for which N q is the surface donor concentration. The hole concentration 
does not appear in Eq. (69) because p n (0) « N q or p R (0) ^ 0. 

2.2.8 Built-in Fields - High and Low Injection Levels 

Built-in fields are usually defined as an electrostatic field 
intensity that exist under thermal equilibrium [64-71]. In treating 
high injection cases, the intensity of these built-in fields may be 
markedly affected [72-76]. Moreover, the high injection level may 
establish another field component that is attributed to photoexcited 
carriers which is strikingly similar in form to the built-in field that 
may be present for immobile charge impurities in thermal equilibrium. 
Therefore, in the following discussion the built-in field label is 
extended to nonequilibrium cases as well. 

The built-in field components will be treated using the ambipolar 
method [72,73] which may be applied for the entire range of injection 
levels (i.e., low to high). This approach is the most direct method to 
illustrate high injection effects on the fields present under equilibrium 
and nonequilibrium conditions. 
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The electric field intensity that may be present in the p-type 
base region is given by 


where 


E pj = E paj + E pBNj + E plPj + E pPEj 


b p.i J l 


■P2J q^jPpj 


E = (b - E )_£2i F 
pBNj ^ pj ^P r pj L pBNoj 

F - kT 2 dn iej 
pBNoj q n iej . dx 


: pIPj 


n n 

(b . - E ) + E 

pj o' r . ^ p r . 

PJ ■ PJ w PJ 


■pIPoj 


-pIPoj 


a_L®si 

q dx 


, t b . - t dn 
: = KL PJ P Pej 

'pPEj q T p . dx 


r w 5 (b qj + { p>"po + «p n aj 


PJ 


!l/2j~fG£V 

F -l/2 (_e Gj~ n j ) 


b ,.-! - n J 


MJ = 


PJ 


°1 = J n M + V x ) 


(70) 

(71) 

(72) 

(73) 

(74) 

(75) 

(76) 

(77) 

(78) 

(79) 

(80) 


Each of the field components in Eq. (70) and the parameters given in 
Eqs. (71) to (79) are identified by the subscript j, which denotes the 
j mesh point. Each is evaluated at the identified mesh point. In Eq. (70), 
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E . represents the field required to conduct the total current ( ) , 
pnj 1 

EpBNj t ^ ie e ^ ectnc Present in nonequilibrium and E pBNoj . in 

equilibrium attributed to bandgap narrowing, E pIp j is the field present 
in nonequilibrium and j p 0 j i n equilibrium due to the impurity concentration 
profile, and E ppE j is the field established by the distribution of photo- 
excited minority carrier electrons. 

Eq. (70) is obtained from a bipolar carrier treatment, where the 
drift and diffusion components of both electron and hole current contri- 


butions are added. Solving for E pj . and substituting for p p from the 
charge neutrality condition results in Eq. (70). The value of the 


parameter £ p = 1 when the ambipolar treatment is required for high 
injection level, and . £ = 0 for low level injection. 

Eqs. (72), (74), and (76) contain r . in their denominators. 

r J 


Under high injection conditions, n.(= n . + n.) » N.. 

pj poj pej hj 


, the field 


components E pBN j and E p jpj decrease with increasing n p ^ , and when 
n P oj/r p j 5 N Aj /r pj << 1 5 they are eligible compared to E pQ . and E pREj . . 

It is clear that the field components, E pB ^ o j and E pIpoj . , that are present 


in thermal equilibrium do not play a major role under high injection 


levels. With increasing injection level, the ohmic field component, E . , 

PWJ 

ultimately approaches a constant value and E ppEj . approaches an asymptotical 

value because r . -* b .n .. 

PJ uJ pej 

Phenomenologically, the reduction in the role of E pB ^. and 
E p ipj may be understood by associated high injection conditions in a 
semiconductor with that of metallic-like behavior. Under these conditions, 


a built-in field established by the distribution of impurities and their 
effect on the bandgap affects the slope of the conduction and/or valence 
bandedges. There is a tendency to neutralize these field components 
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because high current values cannot be supported. Therefore, the high 
concentration of electrons and holes neutralize these field components. 

On the other hand, in order to support the high terminal currents under 
hi qh solar concentration, an electric field, E ^ . , of moderate value 
must be present. The corresponding voltage that is required to conduct 
this current is subtracted from the generated photovoltage across the 
depletion region (x^, x^). Thus, this field is self-limiting and the 
terminal current cannot increase without limit. Similarly, the electric- 
field, Epp£j> is also self-limiting, approaching an asymptotic value. 
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3.0 SUMMARY 


The general method to solve the transport equations is discussed. 
A number of the more important phenomena submodels are presented, and 
the reasons for selecting the particular form are discussed. 
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